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Abstract 

In this paper we investigate the possibility of constraining the hypothesis of a fifth force at the length scale of two 
Earth's radii by investigating the effects of a Yukawa gravitational potential on the orbits of the laser-ranged 
LAGEOS and LAGEOS II satellites. The existing constraints on the Yukawa coupling a, obtained by fitting 
the LAGEOS orbit, are of the order of \a\ < 10~ 5 — 10~ 8 for distances of the order of 10 9 cm. Here we show 
that with a suitable combination of LAGEOS and LAGEOS II data it should be possible to constrain a at a 
level of 4 x 10~ 12 or less. Various sources of systematic errors are accounted for, as well. Their total impact 
amounts to 1 x 1CP 11 during an observational time span of 5 years. In the near future, when the new data on 
the terrestrial gravitational field will be available from the CHAMP and GRACE missions, these limits will be 
further improved. The use of the proposed LARES laser-ranged satellite would yield an experimental accuracy 
in constraining a of the order of 1 x 10 -12 . 
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1 Introduction 

In this paper we investigate the constraints that can be posed on the existence of a possible fifth 
force by the analysis of the laser-ranged data to LAGEOS II and LAGEOS Earth satellites. 
We will consider a potential energy including a Yukawa term of the form [Ohanian and Ruffini, 
1994; Ciufolini and Wheeler, 1995; Nordvedt, 1998] 



where G is the Newtonian gravitational constant, M is the mass of the central body, m is the 
mass of the orbiting test particle, r is the distance between the two bodies, a = JV? contains 



In checking the nature of the fifth force it is of the utmost importance to perform experiments 
spanning the widest range of length scales as possible [Nordvedt, 1998]: for experiments at 
laboratory scale see [Krause and Fischbach, 2001] and references therein. Using the orbits of 
LAGEOS and LAGEOS II satellites implies that we are testing the hypothesis of the fifth force 
at a length scale of almost two Earth radii, i.e. 10 4 km: at this scale the constraints on a are 
of the order of \a\ < 10 -5 — 10~ 8 (see Fig. 3.2 (a) of [Ciufolini and Wheeler, 1995]) and are 
derived from a data analysis of LAGEOS. 

Our analysis includes also an evaluation of the error budget in order to account for various 
systematic errors induced by several classical aliasing forces. We will show that it is possible to 
improve sensibly the present limits by using suitably the data from the existing LAGEOS and 
LAGEOS II satellites and the present or near future knowledge of the terrestrial gravitational 
field whose uncertainties represent, as we will see later, the main sources of systematical errors. 

The paper is organized as follows. In section 2 we derive the effects of the Yukawa gravita- 
tional potential on the orbit of a test body with the standard technique of the Gauss perturba- 
tive equations for the rates of change of the Keplerian orbital elements. In section 3 we apply 
the results obtained in section 2 to the Earth-LAGEOS system and discuss the constraints 
posed on a by using a suitable combination of the residuals of the perigee of LAGEOS II and 

^dn general, K and k are not proportional to M and m. 




(1) 



the couplings K and k of the new force to the two bodies^ and A is the finite range of the new 



force. 
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the nodes of LAGEOS II and LAGEOS. The effects of various sources of systematical errors 
are investigated. The role of the proposed LARES satellite is considered as well. Section 4 is 
devoted to the conclusions. 



2 The orbital effects of the Yukawa perturbation 



The acceleration felt by a test body orbiting the mass M in the potential energy given by eq. (|l|) 
is 

GM A „ , , / 1 1 \ A 
a = a + a pcrt = — — r - a&M [ — - — ) r. (2) 



A 1 



In obtaining eq.(^) it has been assumed that ~ 1 — as it should be the case for an Earth 
orbiting satellite with r ~ 10 7 m. In eq.(@) f is the unit vector pointing from the central mass 
to the orbiting body. 

The second term of the right-hand side of eq.@ can be considered a small, central pertur- 
bation of the Newtonian monopole acceleration. It may be interesting to note that, since, in 
general, a may vary from a body to another in the field of the same mass 2 , eq.([J) implies that 
different bodies may be accelerated differently in the field of M [Nordvedt, 1998]. 

Let us work out explicitly the effects of the Yukawa perturbing acceleration on the orbit of 
an artificial satellite. We will adopt the standard approach based on the Gauss perturbative 
equations and the projections of the disturbing acceleration R, T, N onto the radial, along- 
track and cross-track mutually orthogonal directions [Milani et al, 1987], respectively. The 
Gauss equations are 
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2 This happens if k is not proportional to m. 
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dM 2 r f- du ,dU\ 

— — = n R- - VI - e 2 — + cost— , 8 

at na a \ at at J 

where a, e, i, Q, uj and M are the satellite's semimajor axis, eccentricity, inclination, longitude 
of the ascending node, argument of perigee and mean anomaly, respectively. Moreover, p = 
a(l — e 2 ), / is the true anomaly and n = V GMa~ 3 is the Keplerian mean motion. As can be 
noticed from eq.(0), the Yukawa disturbing acceleration has only the in-plane, radial component 
R, so that it can be straightforwardly inferred from eqs.([5])-(||) that the out-of-plane Keplerian 
orbital elements like the inclination i and the longitude of the ascending node Q are not affected 
by it. By evaluating the Yukawa acceleration on the unperturbed Keplerian ellipse, for which 

_ a(l-e 2 ) 

r ~l + ecosf () 

inserting it in eqs.(|3])-(|8|) and, then, averaging them over an orbital revolution, the long-period 
rates of change of the Keplerian orbital elements can be obtained. It turns out that, by 
neglecting terms of order 0(e n ), with n > 2, in the satellite's eccentricity, only the perigee 
uj and the mean anomaly j\A are affected by long-term Yukawa perturbations. Indeed, from 
eq.(||) and eqs.(|?])-(||) it can be obtained 

du n , , 

dM n 2GM(l-e 2 ) 

"df = n + a Jl^)- a na\* ■ (U) 

3 The constraint on a 

In order to constraint effectively the Yukawa coupling a, let us focus on eq. (|T0D . It tells us 
that the Yukawa perturbation induces on the perigee of a near Earth satellite a secular rate 
which, for LAGEOS II, is proportional to 3.06658517 x 10 12 milliarcseconds per year (mas/y) 
via the coupling constant a. What is the constraint posed on a by the experimental accuracy 
with which it can be possible to measure the perigee rate? In the case of this Keplerian orbital 
element the observable quantity is r = eau. So, by assuming an experimental error of, say, 
^ r exp = 1 cm over 1 y ear > f° r LAGEOS II, which has e = 0.014 and a = 1.2163 x 10 9 cm, we 
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have tf^exp = 12 masQ. This yields a relative accuracy on a 

4 x KT 12 . (12) 



'exp 

exp 



5a 
a 



It is interesting to note that the same estimate for the proposed LARES satellite [Ciufolini and 
Matzner, 1998], which is planned to have a larger eccentricity, Clares — 0.04, would yield a 
relative accuracy of almost 1 x 10~ 12 . 

However, this estimate does not include any systematic errors. As it is well known from 
various proposed or performed tests of General Relativity with LAGEOS satellites [Ciufolini 
et al, 1997; 1998; Ciufolini, 2000; Iono and Pavlis, 2001; Iono, 2001; Iorw et al, 2002; Iorw, 
2002], in such kind of measurements there are lots of competing classical forces which may 
act as superimposed biases affecting sensibly the precision of the measurements. Then, the 
evaluation of the systematical errors induced by them is of the utmost importance. 

The main source of systematic errors is represented by the mismodelled precessions of the 
even zonal harmonics of the static part of the geopotential. In particular, the first two even 
zonal harmonics Ji and J4 are the most insidious. In order to cancel their impact, as proposed in 
the PPN LAGEOS experiment [Iorio et al, 2002], the following combination of orbital residuals 
could be usedQ 

6dj u + Cl 5n n + ddd 1 = axYnk, (13) 

where 

c x = -0.86, (14) 
c 2 = -2.85, (15) 
x Y uk = 3.06658517 x 10 12 mas/y. (16) 

The coefficients of eq. ([T3|) depend on the orbital parameters of LAGEOS and LAGEOS II and 
are obtained in order to cancel the contributions of the first two even zonal harmonics of the 



3 The other existing passive geodetic laser-ranged satellites are unsuitable because their eccentricities are 
smaller than LAGEOS II, except for Starlette whose orbit, however, is known less accurately than that of 
LAGEOS II for various reasons. 

4 Notice that, since LAGEOS enters the combination of eq.(|l3|) only with its node fl, and since the Yukawa 
perturbation does not affect such Keplerian orbital element, cq.([l3|) is insensitive to the possible differential 
falling of LAGEOS and LAGEOS II in the terrestrial gravitational field. 
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geopotential to the measurement of a. It is intended that the residuals would account for the 
Yukawa perturbation in the sense that it would be viewed as an unmodelled feature not included 
in the force models adopted in fitting the satellites' orbits. According to the covariance matrix 
of the most recent available Earth gravity model EGM96 [Lemoine et al, 1998], the systematic 
error induced by the uncancelled even zonal harmonics of the geopotential amounts toQ 

zonals 

| ~7xl(T 12 . (17) 



5a 
a 



Regarding the time-dependent part of the Earth gravitational field, the estimates of [Iorio et 
al, 2002], adapted to this context, yield for an observational time span of 5 years 

harmonics 



5a 
a 



1 x 10~ 12 . 



The most relevant non-gravitational perturbations are the direct solar radiation pressure and 
the Earth's albedo. Their impact can be evaluated from [Lucchesi, 2001] for 5 years as 

/ 5o\ non— s rav 



a 



x 10~ 12 . (19) 



The errors induced by the direct solar radiation pressure and the Earth's albedo have been 
added quadratically, as suggested in [Lucchesi, 2001]. Regarding other subtle non-gravitational 
perturbations of thermal origin acting on the orbits of the LAGEOS satellites, their effects on 
the perigee of LAGEOS II are currently under accurate evaluation. 

Then, a reliable estimate of the total systematic error induced by various classical pertur- 
bations on the measurement of a through eq.(|13|) over a 5 years time span yields 



' r- \ total 

5a \ 



1 x 10~ 11 . (20) 



\ a J 

In obtaining eq. ( |2"0"D we have summed in a root-sum-square fashion the gravitational and non- 
gravitational errors assumed to be independent. Instead, the gravitational error due to the 
static and time-dependent parts of the Earth's gravitational field have been simply summed 
up in view of a reciprocal correlation. 



5 This estimate has been obtained by considering the geopotential harmonics up to degree I = 20. This is well 
justified by the insensitivity of LAGEOS satellites to the higher degree terms. Moreover, this fact makes our 
estimate reliable because the higher degree terms of geopotential in EGM96 are not particularly well determined. 
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4 Conclusions 

In this paper we have shown that by using a suitable combination of the orbital residuals of the 
perigee of LAGEOS II and the nodes of LAGEOS II and LAGEOS laser-ranged Earth satellites 
it would be possible to constrain effectively the Yukawa coupling a of a possible fifth force at a 
length scale of almost two Earth's radii. The experimental sensitivity in measuring the perigee, 
which would be affected by the Yukawa force, of LAGEOS II would allow a precision on a of 
the order of almost 4 x 1CT 12 over a time span of 1 year. According to the present knowledge 
of the terrestrial gravity field, the constraint posed by the systematic errors is of the order of 
1 x 1CT 11 over an observational time span of 5 years. These estimates should greatly improve 
when the new and more accurate data on the Earth gravity field from the CHAMP and GRACE 
missions will be available in the near future. The use of the proposed LARES satellite, with 
its larger eccentricity, would allow to improve the experimental constraint to 1 x 1CT 12 . 
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